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Abstract 

Using a microscopic numerical approach suitable to describe disordered 
antiferromagnets, with application to FexZni-xF2i it is shown that the 
characteristics of the spin glass phase found for a: = 0.25 is much in agree- 
ment with the scenario predicted by the scaling theory of the droplet model 
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The Edwards- Anderson (EA) model [1] has become the prototype model used to 
perform most theoretical studies of spin glasses (SGs). The central feature here is 
to capture the two essential ingredients, namely disorder and frustration. The latter 
is introduced by competing ferromagnetic (F) and antiferromagnetic (AF) interac- 
tions, while the former requirement is built into the model by allowing the exchange 
couplings to have a given probability distribution. A simple solution[2] of the model 
was attempted by assuming infinite-range exchange interactions (SK model), but it 
was later shown[3, 4] to be a quite hard and complex task. This rather unusual mean- 
field theory is characterized by ultrametric-based hierarchy, with infinite degenerate 
ground states unrelated by any obvious (if any) symmetry. This picture proved 
to be of undoubted importance in several fields such as biology and optimization 
problems [5], where large effective coordination numbers mimic a high-dimensional 
space. The question we address, however, is that related to the original motivation: 
is the mean-field theory a good starting point for real short-range (Ising for sim- 
plicity) three-dimensional SGs? This question gave rise to several controversies as a 
more phenomenological approach, namely the "droplet model" [6, 7] , offered a quite 
distinct alternative scenario. The droplet model is based on numerical results of the 
short-range EA model, supplemented by scaling and renormalization-group ideas. 
In particular, it predicts a nondegenerate ground state, except for the trivial up and 
down possibilities. Moreover, in the droplet picture the SO phase is destroyed by 
an applied external magnetic field {H) and thus, contrary to the mean-field pre- 
diction, no equilibrium de Almeida-Thouless (AT) line exists. Experimental [8] and 
theoretical studies [9] facing these controversies have found it difficult to announce 
clear-cut results, in part because of the very slow dynamics plaguing these systems, 
a characteristic shared by both competitive descriptions. 

In this work we study the above-mentioned problems using a direct microscopic ap- 
proach suitable to describe short-range three-dimensional insulating Ising SGs. A 
random-site dilution is thus used instead of the random-bond approach of the EA 
model. In fact, we shall parametrize our model Hamiltonian in the expectation that 
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a description of the experimental observations [10] on the compound Fe^Zni 3.^2, 
as a function of dilution x, is achieved. 

Very recently, we showed[ll] that the presence of a very small frustrated interac- 
tion is the mechanism underlying the appearence of the SG phase at if = in 
-^60.25-2^^0.75-^2- Here we extend our study and calculate the spin-spin correlation 
function in the ordered SG phase and the characteristics of the field-induced do- 
mains both in the SG and in the "random field Ising model" (RFIM) phases. While 
the RFIM phase is very stable under an applied field, its effect on the SG phase 
is dramatic regardless of the value of the field. Our results are in agreement with 
the scenario predicted by the droplet model, including the chaotic behavior of the 
correlation function as the spin distance varies at constant T. The main quanti- 
tative result is the value of the exponent governing the H = SG ordered phase, 
yT = 0.19 ± 0.03. Moreover, we show that an irreversibility (AT) line exists re- 
gardless of the presence of frustration. We conclude that these lines are metastable 
lines, as also supported by the droplet model, and that the SG phase is indeed de- 
stroyed by an applied field, as recently experimentally reported for the compound 
Feo.sMno.sTiOs [12]. 

Formal procedure: we perform numerical simulations of the iterative set of equa- 
tions of the local (site-by-site) mean fields (LMFs) derived from the microscopic 
Hamiltonian 

= Yl Je^i^i+SeSiSi+Se ~ f^oY^iSiH, (1) 

<i,5i> i 

where Si = ±2 (as in -Fe"^^), £^^ = 0,1 and i is summed over the three nearest- 
neighbor exchange interactions Ji of the centred tetragonal lattice of FexZni-xF2 
[10]. We choose the values of the exchange constants by keeping unaltered the 
experimental ratios ji = J1/J2 = —0.013, J3 = J3/J2 = +0.053 and such as to fix 
the Neel temperature of the pure system, Tj^(x = 1) = 77. 8K. This LMF technique 
has proven [13] very efficient in describing inhomogeneous and disordered systems. 
The thermally averaged local spin m^, results from the free energy minimization in 
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the form {£i = 1): 

mi=< Si>T=2tQ,n\i{2hi/KT}, (2) 
where the local field at site i is given by 

= ~o Z] Ji+Se^i+5e''^i+6t + f^oH. (3) 
^ Si 

Introducing Cij(T) =< SiSj >t —mirrij^ we find the LMF expression 

<SS y- + ^■^•)/-^^] - cosh[2(/., - h^)/KT] 

' ^ cosh[2(/ii + hj) /KT] + cosh[2(/ii - /i^) /iCT] ' ^ ' 

The spin-spin correlation function of interest is then Gririj) = 
< Cf^ >av , where av indicates an average over distinct random spin configurations, 
such that < Cij >av= 0. The staggered magnetization is defined by 

M5(^,a;) = (2/7V)E^^m^, (5) 

i 

with i summed over a given sublattice, and the SG-EA order parameter by 

QiT,x) = il/N)Y^ej<S]>T^ (6) 
j 

with j taken through the whole lattice. 

We choose a random initial configuration {Si, Si} in the high-T paramagnetic phase 
and average over 50 independent runs. The system is cooled by an amount AT ( 
we take AT = O.ObK for the (H,T) phase diagram , and AT = IK otherwise ), 
using (2) and a proper convergence criterion[ll, 13], down to T = 2K and heated 
back using the same ammount of AT. We take N = 2 • 30^ sites and use periodic 
boundary conditions. 

Numerical data and analysis: first we point out that our data correspond to a proper 
time scale in which a quasi-equilibrium state is achieved ai H = 0. In Fig. 1(a) we 
illustrate a SG orderered phase configuration through part of a randomly chosen 
sheet: x = 0.25, T/Tm{x = 1) = 0.18 and H = 0. The local-field distribution corre- 
sponding to this state is shown in Fig. 1(6). As evidenced from the inset, competitive 
local fields exist only in the presence of frustration. In fact, the effect of frustration 



manifests itself only in the high dilution regime ( a: < 0.5 ) and causes domains of 
reversed spins surrounded by domain walls already when the system behaves as a 
RFIM ( X = 0.48 ) [11]. This spin reversal effect completely changes the balance 
of forces, thus causing the appearence of competitive local fields, particurlarly near 
the percolation threshold in which case a full SG phase emerges. 
In the scaling theory[6, 7] of the droplet model the SG ordered phase is character- 
ized by the exponent associated with the thermal eigenvalue , {T{L)/ J{L)) ~ L~^^, 
of the T = fixed point, where J{L) is an eflFective exchange interaction. It also 
measures the scale of the free energy, F = Y(T)Ly'^, associated with the thermally 
activated droplets of reversed spins surrounded by domain walls. The average size 
of the domains increases as the temperature decreases, with Y{T = 0) ~ J. We 
have measured the correlation function G^, for x = 0.25, and assumed a scaling de- 
pendence, GriT) = f{r) exp(— r/^), where ^ ~ T^^/vt)^ valid in the low-temperature 
region of the H = SG ordered phase. The fitting of the data is shown in Fig. 2: 
we find ut = 0.19 ± 0.03, in very good agreement with the estimates[6] using T = 
renormalization-group treatments of the short-range EA model. The data are for 
pairs of spins ( ten pairs in this case) inside the same domain. For pairs with spins 
in different domains the correlation is null. Notice that the T-behavior of Gr{T), 
shown in the inset of Fig. 2, confirms a finite-T phase transition. In Fig. 3(a) we plot 
the temperature dependence of < SiSj >t , with sites i and j randomly chosen on 
the lattice. Different behaviors are seen, including a change of sign, indicating spin 
reversal by thermal activation of domains. On the other hand, the chaotic behavior 
of < SiSj >T as a function of spin distance at T/Tn{x = 1) = 0.18 is nicely displayed 
in Fig. 3(6). The Lyapunov exponent associated with this bevavior is predicted[7] to 
be = ds/2 — yT, where ds is the fractal dimension of the domain walls. ( is exactly 
one in d = 1 and expected to change very little with the system dimensionality. 
Last, we should emphasize that below T/Tn{x = 1) = 0.18 there are practically 
no spin reversal and thus for all practical purposes the spin configuration of Fig. 
1(a) is the H = SG ground state for a given initial configuration {Si, Si} in the 
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paramagnetic phase; {—Si, Si} leads to a ground-state configuration with all spins 
reversed. Moreover, for a quenched disorder configuration {Si}, the ground state 
remains the same regardless of the spin paramagnetic configuration {Si}, i.e., the 
SG state is self-averaged[9]. 

We now turn to study the system in the presence of an applied fied. In a field-cooled 
(FC) cycle random fields of strong magnitudes induce SG-like states, regardless of 
the presence of frustration. These spin structures are formed by domains of the two 
(up and down) equivalent H = SG-ground states[6, 7], i.e., the original configura- 
tion of Fig. 1(a) breaks up into domains of length scale, Lh ~ ( J/if)^/'-'^/^^^^'*, which 
independently align with the field. It means that the field is a relevant scaling vari- 
able governed by the (positive) exponent yn = d/2 — yx, {H{L)/ J(L)) ~ L^". The 
spin structures of these SG-like states are illustrated in Fig. 4(a) by indicating the 
domains of reversed spins (□) relative to the H = SG ground-state (Fig. 1(a)), for 
X = 0.25, T/Tisf{x = 1) = 0.18 and H = LOG. The associated distribution of local 
fields are shown in Fig. 4(b). The inset of Fig. 4(b) shows that, in the absence of 
frustration, a zero-field-cooled (ZFC) cycle induces a metastable long-range highly 
nonuniform AF order as a result of the local AF interactions predominating in the 
system. We should stress that we have access to the microscopic local spin values 
on our chosen sheets. It is thus observed that the local spin values on the structure 
of Fig. 4(a) and on the mentioned AF structure (inset of Fig. 4(b)) are the same of 
those in Fig. 1(a) ( for any H < lO^G ), except for a change of sign for those spins 
suffering a spin reversal and a small magnetization component in the field direction 
( experimentally called[10] Mpc and Mzfc, respectively). This is so because our 
fields states were "prepared" using the same initial configuration of Fig. 1 (a) . 
The characteristics of these field-induced states can also be studied by measuring 
the effect of frustration on the EA-SG order parameter. In Fig. 5 (a) we see the 
qualitative change of behavior as a result of an applied field in a FG cycle: while in 
the H = state frustration is a fundamental ingredient ( Q = if = 0), under 
a field, regardless of its value, frustration decreases the Q value ! Saturation of Q 
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is achieved at about lO^T. To further evidence the special character of the H = 
SG phase found for x = 0.25, we have also studied the effect of frustration on the 
system for x = 0.48, a regime dominated by random-field effects on top of the long- 
range AF ordering. In Fig. 5(6) we see that frustration always decreases Q, even for 
H = 0. The effect of frustration peaks at an intermediate temperature, because in 
this situation neither full ordering nor full disorder predominates in the system. 
Finally, in Fig. 6 we display the field-induced AT - irreversibility lines (•) obtained 
by the loci of the points in which Mpc differs from Mzfc- We see that exponents 
compatible with the experimental results are found regardless of the presence of 
frustration. As we see from the insets of Fig. 6(a) and 6(b) the FC (□) and ZFC (o) 
correlation functions, for H = 0.5G, peak at a temperature below the irreversibility 
lines, much in agreement with the H = correlation function in the presence of 
frustration. The data are for sites in the same domain. In the absence of frustra- 
tion, however, the H = correlation function is null (inset of Fig. 6(a)). These 
results show that in our model for a highly disordered magnetic system, metastable 
"glassy" states, characterized by the domains of up and down H = spin configura- 
tions of a SG state at a fixed temperature, are induced by an applied field as a result 
of the strong random fields generated in the system [14]. The effect of frustration is 
only relevant at very low fields, in which case the approach to SG criticality is only 
felt by the frustrated system. Though a specific low-temperature SG configuration 
depends on a given initial paramagnetic configuration {Si,8i}, the ordered phase is 
uni vocally characterized by the thermal eigenvalue yx of the T = fixed point. 
In conclusion, we have presented results based on a local mean-field numerical simu- 
lation of disordered AFs, with application to FexZni-xF2, describing several features 
of these systems as a function of concentration. In particular, for a concentration 
near the percolation threshold, the system exhibits a SG phase with a scenario in 
agreement with the scaling theory of the "droplet" model. 
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Figure Captions 



Fig.l. (a) Part of a SG "ground-state" spin configuration: x = 0.25, T/Tn{x = 1) = 
0.18 and if = 0; O (n) indicate up (down) spins; o indicate zero local spins, 
(b) Local field distribution of the SG "ground state" . 

Fig.2. Best fitting of the low T-dependence of the correlation function: Gr{T) = 
/(r)exp(-r/0 , <^ ~ T^^vt) ^ yelding yr = 0.19 ± 0.03 . The inset shows the 
peak of Gr{T) at the freezing temperature Tf. 

Fig. 3. (a) T-dependence of < SiSj >t , with sites i and j randomly chosen on the 
lattice, (b) Ghaotic behavior of < SiSj >t as a function of spin distance at 
T/Tn{x = 1) = 0.18 . 

Fig.4. (a) Part of a FG spin configuration: x = 0.25, T/Tn{x = 1) = 0.18 and 
H = IG; n indicate spin reversal relative to the SG "ground state" of Fig. 1(a) 
. (b) Local field distribution of a FG spin configuration. The inset shows a 
ZFG local field distribution in the absence of frustration. 

Fig. 5. Effect of frustration on the T-dependence of the EA order parameter for several 
values of fields (magnitude in Tesla): (a) FG cycle for x = 0.25; (b) ZFG cycle 
for X = 0.48 . 

Fig. 6. Field-induced AT irreversibility lines: dashed lines indicate experimental es- 
timate of the crossover exponent, ^ 3.4; dot-dashed lines indicate the AT 
value, = 3.0; full lines are best-fitings. The insets show the T-dependence of 
the correlation function: • indicate data at if = 0; □ (o) are FG (ZFG) data at 
= 0.5T; (a) in the absence of frustation, = 3.6 ± 0.9 ( notice that GriT) is 
null at if = ); (b) in the presence of frustration, (j) = 3.8 ± 0.6. 
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